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It has been suggested recently that in the presence of suitably “warped” extra dimen¬ 
sions, the low-energy limit of pure gauge held theory may contain massive elementary 
vector bosons localised on a “brane”, but no elementary Higgs scalars. We provide 
non-perturbative evidence in favour of this conjecture through numerical lattice mea¬ 
surements of the static quark-antiquark force of pure SU(2) gauge theory in three 
dimensions, of which one is warped. We consider also warpings leading to massless lo¬ 
calised vector bosons, and again End evidence supporting the perturbative prediction, 
even though the gauge coupling diverges far from the brane in this case. 


CERN-TH/2002-319 
January 2003 



1. Introduction 


In standard Kaluza-Klein dimensional reduction of pure gauge theory, the original (say, 
five-dimensional) theory has an effective description in terms of a four-dimensional 
theory, whose lightest degrees of freedom are in the Coulomb or conhnement phase 
(depending on the group), and have a wave function spread out evenly in the fifth 
dimension. It has recently been demonstrated PP that if the fifth dimension is suitably 
“warped”, this pattern could change qualitatively, at least in the Abelian case: the 
low-energy dynamics can still be four-dimensional, but now with massive elementary 
vector bosons, and with a localised wave function along the extra dimension. Thus, 
extra dimensions could potentially provide an alternative for the Higgs mechanism. 

While there is no doubt about the viability of this mechanism in the Abelian case, 
where all computations can be carried out analytically, things are more complicated in 
a non-Abelian theory. For a specific choice of the warp factor the low-energy effective 
action looks much like a four-dimensional gauge theory, but with a gauge non-invariant 
mass term. For an Abelian case, this is still a renormalisable theory, whereas for non- 
Abelian groups it is in general not (see, e.g., ref. j2]). This means that the heavier 
modes cannot decouple from the low-energy dynamics. The hope is that they might 
nevertheless only introduce small contributions, like higher order operators do in chiral 
perturbation theory, but this has so far not been demonstrated explicitly. 

Another way to express the problem is that a gauge theory with a mass term “in¬ 
troduced by hand” may be considered the inhnite Higgs-mass limit of a gauge-Higgs 
theory with spontaneous symmetry breaking la-a, and is therefore strongly coupled, 
at energies of the order of the vector boson mass. Therefore the viability of perturbation 
theory must again be checked by non-perturbative means. 

There are also other types of warp factors, discussed in connection with the localisa¬ 
tion of gravity [5] and gauge fields on a brane, which lead again to a lower dimensional 
effective theory, but this time with massless vector bosons (see, e.g., )• This 

requires asymptotically small warp factors (or, in other terms, large gauge couplings) 
far from the brane As in the previous case, the validity of perturbation theory is 
then in question. Some aspects related to this mechanism were already studied with 
numerical methods in in|. 

The purpose of the present paper is to study the issue of strong coupling with lattice 
Monte Carlo simulations. To simplify the analysis, we would like to separate the 
problem of non-renormalisability of the higher dimensional original gauge theory from 
the problem of a large coupling constant far from a brane. To this end one can study 
a compactification from four dimensions (4d) to three dimensions (3d), or even three 
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dimensions to two dimensions (2d). For the practical reasons that less computer time 
is required, and some exact results are available in 2d physics, we choose here the 
latter case. Nevertheless, we should expect the main features to carry over to higher 
dimensions, as well. 

The outline of the paper is the following. We review some basic aspects of the 
mechanism in Sec. El We introduce our observables and determine their behaviour in 
the Abelian case in Sec. El The lattice formulation is presented in Sec.|31 and numerical 
results for the Abelian and non-Abelian cases, in Sec. El We conclude in Sec. El Some 
technical details are discussed in the Appendices. 


2. The mechanism in review 


We start by reviewing the basic properties of the mechanism introduced in P], in the 
Abelian case. The Euclidean continuum action is 

= Jd^^x JdzA{z)^F^,Fi,, , ( 2 . 1 ) 

where Ff^^ = df^A^, — dj^Afx, and fl = 1,..., d -|- 1, ^ = Xd+i- An index without a tilde 
runs as fl = The function A(; 2 ) > 0 breaks the {d + l)-dimensional Lorentz 

invariance. We however assume the special breaking pattern that terms containing 
F^z are multiplied by the same function. We also take A{z) to be an even function of 
z, and refer sometimes to the plane z = 0 as the “brane”. 

We choose units such that Eq. (EH) should roughly correspond to an effective d- 
dimensional action of the form 

~ , ( 2 . 2 ) 

where g is the gauge coupling (which of course plays no dynamical role in the non¬ 
interacting Abelian case). Thus, 


[2lJ = GeV, [F^J = GeV^ |/] = GeV*-", [/a(j)] = 


(2.3) 


To proceed, we assume that one can make the gauge choice A^ = 0, without intro¬ 
ducing any singularities. It should be noted, however, that this may not always be 
the case in a strict sense, in a non-Abelian theory. If for instance the extent of the z- 
direction and A( 2 ;) are finite, like at hnite temperatures, then A^ behaves effectively 
like a dynamical adjoint-charged scalar field related to the global symmetries of the 
system, which can even get spontaneously broken [THl IIHI120] (for a recent study in the 
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context of extra dimensions, see IZQ). For the purposes of this Section, though, this 
possibility can be ignored. It should perhaps be stressed that all the observables to be 
introduced later on, as well as the main lattice simulations carried out, are explicitly 
gauge invariant, so that our conclusions are by construction based on data which is 
independent of the gauge choice. 

We now carry out a mode decomposition of the functional dependence of the helds 


on the ^(-coordinate, 

n 

(2.4) 

Units are chosen such that 


|V-„] = Is] = 

(2.5) 


The real functions ipni^) are assumed to satisfy the second order Sturm-Liouville linear 
differential equation. 
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( 2 . 6 ) 


Here are real, because the differential operator is Hermitean. They turn out also 
to be non-negative. We denote the mode constant in 2 ; (whether normalisable or not in 
inhnite volume) by the index n = c, while general normalisable states with non-negative 
masses are labeled by n > 0, with even (odd) indices denoting states symmetric (anti¬ 
symmetric) in z —>• —z. Explicit solutions of Eq. ()2.6jl for various A(z) are discussed 
in [Appendix A[ Note that if the constant mode is normalisable in inhnite volume, then 
the indices n = 0 and c refer to one and the same mode. 

Together with the normalisation condition 


J A{z)^/Jmiz)^/Jniz) = Smn, (2.7) 

Eq. ()2.6j] guarantees that 

/ = ml5ran- ( 2 . 8 ) 

J z 

Note also that the completeness relation can be written as 


^■iljn{z)ilJn{z') = A ^{z)5{z - z'). (2.9) 


The quadratic part of the action then becomes 

SWD _ j d''x + \miAlA";). (2.10) 
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If the lowest mass is zero or much smaller than the masses of higher modes, we have 
an effective d-dimensional field theory at low energies: it is described by the term with 
n = 0 (or n = c) in Eq. ()2.1()j) . 

Now, if the extent of the z-direction is finite and A(z) is regular, or if 

/ OO 

dz A{z) < OO, (2-11) 

-OO 

then Eq. (ESI) clearly has a normalisable zero mode solution, with ipoiz) 
and rriQ = ml = 0. The normalised form of this solution is 


Then the low-energy effective theory is simply a standard pure gauge theory. The 
condition Eq. (ITTTl) implies that hm^^ooA(^) = 0 and, therefore, that the effective 
d+1 dimensional gauge coupling A“^( 2 ) is large far from the brane. In the non-Abelian 
case this fact may, in principle, invalidate the perturbative arguments just presented, 
and thus provides a motivation for a lattice study. 

If the condition in Eq. (ITTTD is not satisfied, then the constant mode V’c effectively 
decouples (since ^ 0 ); 7 ^ 0 ; and we have massive vector bosons without any 

scalar particles. Furthermore, a mass hierarchy tUq -C mf can be achieved with some 
choices of warp factors (see | Appendix A| and ref. P), provided that A(0)/A(2;o) 3> 1, 
where zo is a point where A( 2 ;) reaches its minimum value. Thus, a large mass ratio 
again only appears if the effective higher-dimensional gauge coupling A“^( 2 ;) is large, 
but now at a hnite distance Zq from the brane. 

Another subtle point with the case mo 7 ^ 0 is that the low energy action is seemingly 
not gauge invariant (see P for a discussion of gauge transformations). In the Abelian 
case the theory is nevertheless renormalisable, even if some interactions were added 
(see, e.g., |23)- This is no longer true for non-Abelian theories, and the question 
appears whether the higher lying modes decouple or not. 


= -ipc constant 

( 2 . 12 ) 


3. Static force in the continuum 

In order to distinguish the two different regimes (with and without the massless vector 
mode 'ipc) we shall employ the standard order parameter for conhnement, the static 
force between two heavy test charges in the fundamental representation. We measure 
the force at a fixed z] for actual mechanisms for the localisation of scalars and fermions 
in the vicinity oi z = 0 see, e.g., jHj and references therein. 
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For now, we shall restrict to d = 2, the case we have actually studied with lattice 
simulations. We consider a rectangular area with (Axi,Aa; 2 ) = and define a 

Wilson loop around the rectangle, 

W{r,t;z) = (RoTtV exp{i z) 

= (ReTrVexp{iJ2i’n{z) ([ A^{x)dXf,)y (3.1) 

-rj 


The static potential can then be obtained as usual, 

V (r; z) = — lim - In IF(r, t; z). 

t^oo t 

A lowest order computation gives 


V{r;z) 




dp - 1 ^ jjljz) / _ 
271 p‘^ + ml n 2m„ ^ 


(3.2) 


(3.3) 


The static potential, itself, is of course not a physical observable. Depending on the 
spectrum m„, its absolute value can be ultraviolet divergent, and in any case sensitive 
to ultraviolet physics. Therefore we rather address its derivative, the force F{r-,z), 


F{r; z) 


dV(r; z) 
dr 


(3.4) 


According to Eq. ea, 

= (3-5) 

n ^ 

We note from Eq. (E3) that an external source couples to the mode n via = V’n(^)- 
The signatures expected from F{r]z) can thus be summarised as follows. In the 
case that the zero mode exists, mo = me = 0, the force should approach a constant at 
large r, 

F{r;z)^^ijl, (3.6) 

because massive modes give contributions screened at distances r>l/mn. On the 
other hand, in the case of interest to us where mo ^ 0, mo <S mi, and the zero mode 
decouples (-^c —>■ 0), we expect 


F(r; z) 



(3.7) 


It is thus our objective to show that the force does get screened, but only on large 
distances, as determined by 1/mo. 
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While we focus on the force in this paper, a behaviour qualitatively very similar can, 
particularly in the Abelian case, be obtained from various correlators of local gauge 
invariant operators. For completeness, we discuss one example in [Appendix B 

So far we have discussed the potential in the Abelian theory. In the non-Abelian 
case, the self- and cross-interactions between modes make obviously a fully analytic 
computation impossible. However, if dimensional reduction indeed takes place then, as 
discussed in [Appendix C[ the only change in the long-distance force is a colour factor, 
the quadratic Casimir of the fundamental representation, Ca = — l)/{2Nc)'. 

F{r-,z) ^ CAF{r-,z) . (3.8) 

This simple relation, which allows us to directly compare the asymptotic non-Abelian 
force with the Abelian one, is obviously specihc to 2d physics. 

In the non-Abelian case, it is useful to also dehne couplings characterising the cubic 
and quartic self-interactions of the fundamental mode. Let us introduce 

93= [ 9l= [ , (3.9) 

J Z J Z 

and construct the dimensionless quantities 

2 

_ Os _ 5^4 /Q -I 

In order for the low-energy effective theory to be “close” to a gauge theory, these 
numbers had better be close to unity. In particular, if the zero mode is normalisable 
and therefore 'ipo = iIjc-, we have exactly = 04 = 1. In the opposite case of uiq ^ 0, 
we have ^ 1 and 04 7 ^ 1. Thus, the breaking of gauge invariance in the low-energy 
sector manifests itself both through an effective mass term in Eq. (ITTUl) . and through 
non-universal self-interactions which differ from the coupling of the modes to external 
sources. 


4. Static force on the lattice 

As mentioned above, in the non-Abelian case the heavier modes cannot decouple, 
because they are needed to guarantee renormalisability. It is therefore not obvious how 
well the analytical estimates presented in Sec. El really hold. We will hence study that 
system with simple numerical lattice Monte Carlo simulations. 

In fact, to account properly for finite size and finite lattice spacing effects, we will 
carry out small scale simulations for the Abelian system, as well. Thus, we can directly 
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compare the two sets of data, with similar volumes and lattice spacings. This may be 
useful because the “sharp” and “smooth” weight functions to be introduced contain a 
small scale hierarchy, which tends to lead either to hnite size or finite lattice spacing 
effects in lattice simulations. Still, both sets of results turn out in most cases to remain 
close to the analytic continuum estimates. 

In the Abelian case, we discretise the action in Eq. m by using the so-called 
non-compact formulation: 

= m) 

^ X jKu ^ 

where a^£;(x) = a^{x) -A- a^{x -|- /i) — + v) — aa^x), ap,{x) = aAjx[x), and a is 

the lattice spacing. For future reference, we also dehne the link matrix, Uf)_{x) = 
exp[zQ;^(a;)]. The dimensionless coupling constant appearing in Eq. (14.111 is taken to be 

I3c(z} = ^ . (4.2) 

a 

In the non-Abelian case, we employ the standard Wilson action, 

4"“' = T,0c(z) E E(l - A (4.3) 

2 X ■‘'C 


where the naive continuum limit implies 




(4.4) 


Rather than j^dz), we will often equivalently refer to Ao/a to fix the lattice spacing, 
where Aq = A{z = 0). Note that we can view a as being constant throughout the 
lattice: in our case a non-constant I3g{z) does not imply varying lattice spacing. 

It should be noted that, as we have discussed in |Appendix A the value of Aq does not 
affect at all the spectrum obtained in the non-interacting limit. For a weak coupling, 
the criteria for discretisation and finite volume effects to be small are simply 


a -C — -C T, T, 
mo 


(4.5) 


where L, T are the linear extents of the system in the r and t directions, respectively. 
On the lattice, however, Aq determines the strength of interactions. In general, lattice 
discretisation effects are larger and the gauge theory more strongly coupled where Pciz) 
is smaller, if arrio is kept hxed. We return to this issue presently. 
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It is useful to note that if we think in terms of the mode decomposed action 
in Eq. (ITTUl) . then the mode n can effectively can be assigned a 2 d action at any 
hxed z, with 




2 iVe 


2iVe 


(4.6) 


Parameterising the dimensionless 2d link matrix as 



_ ^iaipniz)T^Al{x) 


(4.7) 


where are the Hermitean generators of SU(A'"c), the naive discretisation of the n = 0 
part of Eq. (ITTUl) then becomes 


srh)= i: [i: (i -Re Ti' d?) + i: (i - Ti' 


(4.8) 


An action of the form in Eq. (USD is of course not gauge invariant, and thus in general 
not (perturbatively) renormalisable. It also does not yield the correct naive continuum 
limits for the cubic and quartic self-interactions, if eta, 04 7 ^ 1. Nevertheless, we might 
still hope Eq. (USD to contain some qualitative features of the effective low-energy 
dynamics, to the extent that the theory is weakly coupled, and the results are only 
moderately dependent of the lattice spacing (or ultraviolet physics), as may indeed be 
expected to be the case in two dimensions [2]. 

The observable we measure on the lattice is the static force. The dehnitions follow 
Eqs. (EH), (EH, EH, only the Wilson line is constructed by multiplying together the 
link matrices around the rectangle, both in the non-Abelian and in the Abelian cases. 
We determine the force then as 


F(r -|- -a; z) 


V{r + a] z) — V (r; z) 
a 


(4.9) 


Note that in the Abelian case the potential is invariant in r —L — r and the force 
then, for a hnite L, takes the form 


F{r;z) = ^ 

n>0 


V^n(^) 

2 


sinhm„(-| — 
sinh^ 



(4.10) 


instead of Eq. (13.5|) . For the non-Abelian case such a periodicity would only arise for 
a force dehned from the correlator of two Polyakov loops (see, e.g, 0 ). 












5. Numerical results 


We now present our numerical results, obtained with standard Monte Carlo simulation 
techniques. The update is a 1:4 mixture of heat-bath 1211123 and over-relaxation ra 
sweeps. In the SU(2) simulations, we use the following statistical noise reduction steps 
in the Wilson loop measurements: 

1. First, we perform link integration for the links in the t-direction, substituting 
each with the appropriate (and exactly calculable) local statistical average link. 

2. Then, we do two smearing I2SI steps for the links along the (r, z)-plane, “fuzzying” 
the r-sides of the Wilson loops by two lattice units in the 2 ;-direction, with rapidly 
decreasing weights. This enhances the coupling to the lowest modes, which are slowly 
varying in 2 ;. 

Both the link integration and the smearing must be performed taking into account 
the varying coupling I3g{z). We perform between 10^ and 4-10^ sweeps and collect the 
data typically in 100 bins. Errors are estimated with a standard jackknife analysis. 

As our goals here are of a qualitative nature only, we should stress that these are 
still very simple small scale simulations. Presumably our numerics could have been 
signihcantly improved for instance by implementing the advanced methods introduced 
in pn] . 

5.1. Gaussian weight function 

We will start with a study of a Gaussian weight function, 

A(^) = Aoexp(-^m^;2^) . (5.1) 

The spectrum following from it is discussed in Appendix ro and goes as = y/nm, 
n = 0,1, 2,... . The zero mode, 'ipc = does have a finite coupling in this case, since 
J^A(z) is hnite. Note that even though the wave function xjjc is constant in 2 ;, this 
mode is said to have been localised, in the sense that centered around z = 0. 

The original theory has two parameters, Aq, m. As both of them are dimensionful 
([Aq] =GeV“^, [m] =GeV in (2-|-l)d), continuum physics only depends on their prod¬ 
uct. Moreover, the mass spectrum, and thus the dynamics of the Abelian theory, are 
completely independent of Aq. Gorrespondingly, it appears that the non-Abelian the¬ 
ory can be made weakly interacting by choosing a large value of m ■ Aq. This argument 
might fail, however, because mo ■ A{z) is exponentially small at large 

The lattice introduces a further dimensionful parameter, the lattice spacing a. It 
should be chosen small enough such that discretisation effects are harmless. Somewhat 
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Figure 1: The force F{r]z) for the Gaussian weight function, Eq. (I5.1|l . at different 
hxed values of z, in the Abelian (left; volume = 48^ x 14) and non-Abelian cases (right; 
volume = 24^ x 14). The perturbative values are also shown. The Abelian case has a 
hnite slope because of the periodicity discussed around Eq. (I4.1()|l . 


arbitrarily, we then £x a, m, Aq such that 

{arnf = 0.1, 4Ao/a = 60.0. 


(5.2) 


Thereby the theory should be weakly coupled (m ■ Aq ~ 4.7), and also close to contin¬ 
uum behaviour {a ■ m k. 0.3). The extent of the lattice in the z-direction is chosen as 
14a~4.4m“^ (cf. Fig. [T)). 

The massless zero mode present in the system should dominate the physics at large 
distances. Its coupling is independent of 2 ; and, according to Eqs. (Eni), (EHD, 




(5.3) 


where Ca = 1 for U(l), 3/4 for SU(2). The other modes give exponentially suppressed 
contributions, according to Eq. (13.51) . 

In fact, we can easily make an exact continuum prediction for the full force F{r]z) 
in the Abelian case. The analytic values of = \/nm and = 0) are given in 
Appendix mu and can be plugged into Eq. ()4.10jl . The prediction is compared with 
numerical data in Fig. ^ (left). Although the data becomes noisy at large r, we can 
conclude that there is agreement within statistical errors, confirming that discretisation 
effects are under control. 
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Figure 2: The force F{r-,z) for the sharp weight function, Eq. (j5.4jl . at z = 0 in the 
Abelian (left) and non-Abelian cases (right). In the Abelian case (where no noise 
reduction techniques were used) the force is anti-symmetric with respect to r/a = 15. 
For comparison we also show the result from a 2d simulation based on Eq. (gSl). 

The numerical result for F{r-,z) in the SU(2) case is shown in Fig. Q (right). It 
indicates that the large distance behaviour is successfully predicted by the perturbative 
analysis. Moreover, the fact that the constant value to which F{r-,z) tends does not 
depend on z conhrms that the constant mode dominates in that regime. We see also 
the divergence of F(r; z) at small distances, as in the U(l) data. Qualitatively, the U(l) 
and SU(2) cases yield very similar results. This is a peculiarity of 2d physics, however; 
had we compactified onto three dimensions, our expectation would be F{r-,z) ~ 1/r 
for U(l), and yet still the conhning constant force for SU(2). 

5.2. Sharp weight function 

We next study a warped “sharp” weight function, 

A(^) = Aoexp(^—M\z\ + . (5.4) 

The corresponding spectrum is discussed in Appendix IA.2I 

In addition to the requirements for the Gaussian case above, leading to a weak 
effective coupling and small discretisation effects, we are now faced with an additional 
constraint, as well as an additional parameter allowing to satisfy it: we want to tune 
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M/m such that the “fundamental mode” the one with the lightest non-zero 

mass mo, is much lighter than the next mode, with mass mi. This tends to make the 
choice of parameters somewhat less transparent. From the point of view of the infinite 
volume setup, the zero mode V’c is an artifact of the simulation, whose effects ought 
also to be kept numerically small. 

In practice, we choose the parameters as 

am = 0.50, aM = 0.75, 4Ao/a = 35.0 . (5.5) 

The spectrum resulting from these parameters is discussed in Appendix lA.21 The mass 
of the “fundamental mode”, mg, is (mo/m)^ = 0.235. In lattice units, therefore, 

,^o/a = (amo)~^ ~ 4.1 . (5.6) 

The first excited state with a finite coupling 'ipn{z) at z = 0, on the other hand, has 
a correlation length ^ 2 / 0 , ~ 1.4. Thus the fundamental mode should indeed dominate 
the infrared physics. For the lattice size used, 30^ x 18, the effective couplings of the 
zero and fundamental mode are 


a'lpc ~ 0.048, a-^o(O) ~ 0.194 . 


(5.7) 


We observe that because of the finite volume, ai/c is not quite zero yet. The contribu¬ 
tions from the different modes to a?F{r] z) follow from Eq. (j3.5j) and, fortunately, the 
effect of a'lpc turns out to be smaller than our error bars. Note also that a dimensionless 
2 d effective coupling can be estimated as 


V’o 0.79 
_ _ 

vrmo vr ’ 

and the parameters as and a 4 defined in Eq. dSHI evaluate to 


(5.8) 


as = 0.804 a4 = 0.749, 


(5.9) 


suggesting that the perturbative picture, as well as a 2d action of the form in Eq. ra, 
should be qualitatively applicable. 

The data is shown in Fig. |21 In the Abelian case, we insert the values of ipniz = 0) 
and mn computed in Appendix A into Eq. (Id.lOj) . to obtain the continuum prediction, 
shown with the dashed-dotted line. The dashed line shows the contribution of the 
fundamental mode alone. The simulation is quite consistent with the exponential 
decay of F{r]z). However, the periodicity in the r-direction makes the extraction of 
the force quite difficult and noisy when r L/2. 
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In the SU(2) case, the value of the fundamental mode prediction at r = 0 sets the 
scale for the asymptotic confining force, were the fundamental mode massless (dashed- 
dotted line). Clearly, the lattice data is consistent rather with a decaying force, or the 
“breaking of the string”, on a distance scale given by ^ comparison, 

we show also the result for the 2d system, Eq. (gSl), using the tree-level values of 
a'0o(O), amo as input. We observe the same qualitative behaviour, although the (2-1-1)- 
dimensional case leads to a stronger force at small distances, due to the exchange of 
higher modes. 

Finally, we reiterate that because of the finite value of aipc in our finite box (cf. 
Eq. ()5.7|l h the 3d force should at very large distances still approach a finite non¬ 
vanishing value, ~ 0.001, which is however beyond our resolution. 


5.3. Smooth weight function 

We end by studying a “smooth” weight function, 

A{z) = . 

The corresponding spectrum is discussed in Appendix IA.3I 
We choose the parameters and lattice spacing such that 

am = 0.2778, aM = 0.3889, 4Ao/a = 35.0. 

According to Appendix the correlation length of the fundamental mode 

^o/a = (amo)“^ ~ 5.2. 

The couplings of the zero and fundamental mode are, at z = 0, 

aipc ~ 0.054, a'^o(O) ~ 0.134. 

These determine the string tension, according to Eq. (jSHD. We also note 
dimensionless coupling related to the fundamental mode is 

-00 0.69 

_ _ 

vrmo vr ’ 

and the parameters and 04 defined in Eq. (ITTUD evaluate to 

Us = 0.705 0^4 = 0.695, 


(5.10) 


(5.11) 
is then 

(5.12) 

(5.13) 
that the 

(5.14) 

(5.15) 


supporting again the qualitative applicability of the perturbative picture, as well as of 
the 2d action in Eq. (gSl)- 
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Figure 3: The force F{r]z) for the smooth weight function, Eq. ()5.1()jl . at z = 0, in 
the Abelian (left) and non-Abelian cases (right). In the Abelian case (where no noise 
reduction techniques were used) the force is anti-symmetric with respect to r/a = 15. 
For comparison we also show the result from a 2d simulation based on Eq. (gSl). 

The lattice simulation is carried out with the volume 30^ x 18. The data is shown 
in Fig. El We proceed as in the “sharp” case to compute the prediction for F{r;z). 
The U(l) data is consistent with the analytic prediction, although large statistical 
errors make this statement rather weak. For the SU(2) case we again observe that 
the behaviour of the 3d simulation agrees at large distances within statistical errors 
with the massive fundamental mode prediction, as well as with results obtained with 
the 2d action in Eq. (USD, using the tree-level values of 0-00(0)) as input. At 
smaller distances, the static force is stronger in the (2-1-1 )-dimensional system, as noted 
previously. Note also again that because of the hnite value of aipc in our hnite box 
(cf. Eq. (I5.13|l i. the 3d force should at very large distances still approach a hnite non¬ 
vanishing value, ~ 0.001. 


6. Conclusions 

We have studied in this paper some physical properties of a pure gauge held theory, 
living in a space where the coupling constant g‘^{z) = varies along one spatial 
direction. Using a mode decomposition and working in the gauge = 0, an ehective 
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lower dimensional action was already derived in pp. However, in the non-Abelian case, 
several complications arise: all Kaluza-Klein like modes are coupled through cubic and 
quartic terms, possible non-perturbative effects in regions where the coupling is large 
make the validity of the perturbative analysis unclear, and the naively truncated action 
for the low-energy sector is non-renormalisable. These difficulties motivated a lattice 
simulation of the non-Abelian theory in (2-1-1) dimensions, as well as, for calibration, 
simulations of the Abelian theory in the same background (in which case our analytic 
predictions are exact in the continuum limit). 

In the case of a Gaussian profile A( 2 ), our numerical data conhrms the presence of 
a massless constant mode, which gives rise to a constant force both in the Abelian and 
the non-Abelian cases, in spite of the coupling becoming strong at large 2 ;. Note that 
even though the mode is constant in z, this mechanism is conventionally called the 
localisation of massless vector bosons, since A( 2 ;)x[the mode] is sharply centered. 

For two different profiles such that f^A{z) = 00 , on the other hand, where the 
massless mode decouples from the theory, the lattice data is consistent with the long 
distance dynamics being dominated by a single massive localised “fundamental” vector 
mode. Because of a large hierarchy between the mass mg of the fundamental mode, and 
those of the higher modes, this regime can set in even at distances somewhat smaller 
than the correlation length of the fundamental mode. 

Thus we conhrm the qualitative picture based on perturbation theory in these 
(2-|-l)-dimensional systems. It would be interesting to extend the study to a (3-1-1)- 
dimensional case, to check whether our conclusions depend on the peculiarities of the 
two-dimensional effective theory. 
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Appendix A. Energy spectra for various weight functions 

In this Appendix, we determine explicitly the spectra for the various weight functions 
appearing in this paper. 

As mentioned in pQ, we can write Eq. (ESI) in another form by introducing Xn{z) = 
Then the eigenvalue equation takes the familiar form 

- Xn + Vo{z)Xn = mlxn, (A.l) 

where 

VUz) = W%)-W'{z), H/(j) = -AW. (A.2) 

One may also introduce 

V,{z) = W\z) + W'{z) ; (A.3) 

its eigenvalues are the same as those of Vo{z), except that one of the two has a normal- 
isable exact zero mode, and the symmetry properties of the two sets of wave functions 
with the same energy are the opposite Q. Thus, denoting and assuming 

EP = 0, we have = Ei^+^\n > 0. 

It is useful to note that if A{z) = Aq exp(/( 2 :)), then 

IV = -!/', v; = i(/y + l/". (A.4) 

Therefore, the eigenvalues are independent of Aq. 


A.l. Gaussian weight function 

We start by considering the Gaussian weight function, Eq. EH), which implies that 

W = ^rn^z, Vo = z"^ — . (A.5) 

The eigenvalue equation, Eq. (EH), is just of the form of a harmonic oscillator, with 
shifted energy levels, and is immediately solved. We obtain 

rrin = m^/n, n > 0. (A.6) 


Note the existence of a normalisable zero energy solution. It appears here for Vo rather 
than Vs, since J^A{z) is finite. 

We also know the wave functions exactly in this case. 


^n{z) = 


111 

2^/2 



(A.7) 
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Sharp profile, c=1.5 


Sharp profile, c=1.5 



Figure 4: Left: The fundamental mode V’o(^) and the prohle A(^) for the sharp wave 
function, Eq. ()5.4|1 . at c = M/m = 1.5, R = 4.5/m, and a specihc a (see TableH)). The 
horizontal axis is for C, = zm. Right: the hrst excited even modes. 


where Hn are the Hermite polynomials. Using the fact that Hn have the parity of their 
index and that 

= 0) = (-1)" Ldl, (A.8) 

nl 

we arrive at the expression 


'l/2n{z = 0 ) 


1 ^ (- 1 )" 1 
V 2/ n\ a/ttAq ’ 


(A.9) 


where the last step is the large n asymptotic behaviour. This means that the higher 
modes are not only more massive, but also more weakly coupled at z = 0. 


A.2. “Sharp” weight function 

The sharp weight function, Eq. Q, leads to 


W = ^(Msign(z) —m^z) . (A.10) 

It is convenient to rescale everything by m: (/ = mz, c = M/m, = m\/m?. Eq. ()A.1|1 
then becomes 


Xn + 


^(C-csign(C))^ + ^ - C(5(C) 


Xn = 


(A.ll) 
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The solutions only depend on c. Note that although we are using the same notation 
as in Eq. (EH), Xn is here treated as a function of C, rather than z. 

Since the Hamiltonian is invariant under parity, one can classify its eigenfunctions as 
symmetric and antisymmetric under C, —>■ —C- We label symmetric and antisymmetric 
states with even and odd indices, respectively. 

Infinite volnme. We start by discussing the inhnite volume case. Here we solve 
the problem by utilising Vg. Note that, apart from the usual relation 
n > 0, there is now the additional relation that the eigenvalues obtained with Vg and 
Vo, coming with wave functions antisymmetric in [z —*• —z), are trivially related by the 
addition of m^, because the antisymmetric state does not “see” the 6 function at the 
origin. This can be expressed as E^'> = E^'> + m^, for p odd. Therefore, the spectrum 
is of the form 


Ef> 


(A.12) 

ew 


(A.13) 

eP 

= Ef >Ef,... . 

(A.14) 


In other words, Vg has a symmetric state with exactly zero energy, and an antisymmetric 
one with an exponentially small energy, e (cf. Eq. f)A.26fl below). Then W has a 
symmetric ground state with the energy e, and a doublet of states with a much higher 
energy, Ep) Ep) S> e. 

The explicit form of the dimensionless equation with Vg becomes 


~ Xn + 


^(C-csign(C))^-^ + c<5(C) 


Introducing the Kummer function 

a ( a{a + 1 ) 


+ ... = E 


which satisfies 





(A.15) 

(a)n C 

(dn n\ ’ 

(A.16) 


(A.17) 


the general solution of Eq. (jA.lhjl reads, for 7 ^ 0 and denoting —>• u, 


X(0 = e' 


(KI-oV 


AM—-, + B (Id - 


^ 3 (Id -c) 


(A.18) 

where A,B are constants. For ^ < 0, we denote y by xl, with constants 

and for C > 0, XRi with Ar, Br. The symmetric wave functions obviously have Al = 

Ar, Bl = Br, the antisymmetric ones A^ = —Ar, B^ = —B^. 

The boundary conditions we have to impose on the coefficients are: 
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(a) Xi?(0) = Xl(0), 

(b) X'i?(0) = Xl( 0) + cxL,i?(0), 

(c) lim^^ooXfl(C) = lim(;^_ooXL(C) = 0 . 


The second comes from integrating both sides of Eq. (lA.lhfl from —6 to 6 . 

In both the symmetric and antisymmetric cases, the third condition imposes 




r(- 


R, 


(A.19) 


where we used that for large (, 


r(a) 


(A.20) 


For symmetric wave functions, the condition (a) is automatically satished. The 
condition (b) yields 

c0'i . 


Br = 


03 + 

where 0'^ = 0'(—f; y), 03 = 0(^; |; y), 03 = 0X^; |; y)- Combining this with 

Eq. ()A.19|1 . one obtains an algebraic equation for the energy levels: 


(A.21) 


u 3 


r(if^)r 1-1/ 3 




—; — 
2 2 2 ^ 3 


1 — u , ,3 — u 5 


- 0 ( 


2 ’ 2’ 2 'J 


where we made use of 


0'(a;&;C) = -^0(0 +1;^ + i;C) • 


= 0 , 

(A.22) 

(A.23) 


Note that since r(—1 //2) has a pole at i/ = 0, i/ = 0 is a solution for any c, as must be 
the case. 

For antisymmetric wave functions, the condition (a) yields 


Bn = ^An , 

C03 


(A.24) 


with the same notation as above. The condition (b) is then automatically satisfied. 
Together with Eq. ()A.19jl . one again obtains an algebraic equation for the energy levels: 


^(“2)0(“2’ 2’ ^ V^cr(^—)0(^—; 2^ y) = 0 • 


(A.25) 
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In the limit c 3> 1, the approximate solution for the lowest energy level is obtained by 
setting 1 / —0 in the argument of the hrst 0 appearing, whereby 0(0; y) = 1; we 
then get 

\J2/'Kce~^. (A.26) 

Some numerical values for v (expressed as are included in Tabled 

Finite volume. Let us now consider the same system, but in a box with periodic 
boundary conditions at ^ = ±R/m = ±R, rather than in inhnite volume: 

Xl{-R) = Xr{R)- (A.27) 


Then the spectrum changes. We will study this system directly in terms of W, rather 
than Vg. The general form of the solution in terms of the Kummer functions remains 
the same, except that the Hamiltonian at ( ^ 0 has changed by a constant. 

We will now need to impose a boundary condition at ( = R. Integrating the equation 
of motion for ip from R — 6 to R + 6, it is easily seen that 0' must be continuous at 
C = R (because A is). This translates into 

x'r{R) + W{R)xr{R) = 0, (A.28) 

where W = W/m. This means that x’s derivative is not continuous at the boundary, 
which is due to the discontinuity of A' at C = R- The other boundary condition for 
the symmetric states is Eq. ()A.28(1 applied at the origin or, equivalently, condition (b) 
above, 

X^(0) + lE(0)xi?(0) = 0 . (A.29) 

Note that Eqs. (IA.28|1 . (IA.29|1 imply that '0)j(O) = 'iP'r{R) = 0. In the antisymmetric 
case, the complete boundary conditions are 


Xi?(0) — 0, Xr{R) — 0, 


(A.30) 


and Eq. (IA.28jl then imposes that Xr{R) vanish as well. 
For symmetric wave functions, Eq. ()A.29fl implies 


Ar = DB 


R^ 


D = 


(1 - c2)03 + 


C(0'l - 0 i) 

where the notation is as above. Inserting into Eq. (IA.28|) . 

'l - (fl - cf]4>§ +(R- c)Vf + (R- C)(^f - = 0, 


(A.31) 


(A.32) 
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n 

sharp prohle, c 

R = A.b 
lyA) ai/jn{0) 

= 1.5 

R = oo 
„(») 

sharp prohle, c 

.R = 4.5 

lyA) a^/Jn{0) 

= 2.5 

R = oo 
„(») 

c 

0.000 

0.048 

— 

0.000 

0.201 

— 

0 

0.235 

0.194 

0.209 

0.154 

0.148 

0.038 

1 

1.235 

0.000 

1.209 

1.154 

0.000 

1.038 

2 

1.945 

0.116 

1.643 

1.815 

0.101 

1.191 

3 

2.945 

0.000 

2.643 

2.815 

0.000 

2.191 

4 

4.432 

0.113 

3.201 

4.298 

0.106 

2.478 

5 

5.432 

0.000 

4.201 

5.298 

0.000 

3.478 

6 

7.854 

0.112 

4.829 

7.723 

0.109 

3.867 

7 

8.854 

0.000 

5.829 

8.723 

0.000 

4.867 

8 

12.24 

0.112 

6.504 

12.11 

0.110 

5.327 

9 

13.24 

0.000 

7.504 

13.11 

0.000 

6.327 


Table 1: Eigenvalues for the sharp prohle. For the wave functions at origin one needs 
also the values of am, Ao/a, entering as a'0„(O) oc (am)^/^(Ao/a)“^/^, if c, R are kept 
hxed; we have assumed am = 0.5, Ag/a = 8.75, for both values of c. 


where ^ ^ ^ ^ 

|i )- This equation determines n, after which the actual energy level is 
found by adding 1. 

For antisymmetric wave functions, Eqs. (jA.Onjl imply that 


01 


R 

^ Aj,. 


Bjl — ;— Ajl, Bjl — — 

C03 (c - R)(j)f 


(A.33) 


These are compatible only if 


(C - = C030f , 


lR 


(A.34) 


which determines the eigenvalues. 

As an example, some numerical values are given in Table d The doublet structure 
as discussed above is visible in the inhnite volume results (levels 1 & 2; 3 & 4; ...). 
It can also be observed that for small n, the spectrum is roughly linear in n (like in 
Eq. ()A.6fl for m^), while for large values within a hnite box, it starts to resemble more 
the corresponding spectrum ~ n^. 
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Smooth profile, c=1.4 





Smooth profile, c=1.4 





Figure 5: Left: The fundamental mode ipoiz) and the prohle A(z) for the smooth wave 
function, Eq. dd.lOj] . at c = M/m = 1.4, R = 2.5/m, and a specihc a (see Table El). 
The horizontal axis is for C, = zm. Right: the hrst excited even modes. 


A.3. “Smooth” weight function 

The smooth weight function, Eq. (ICTHl . leads to 


W = 



(A.35) 


We rescale again everything by m: ( = mz, c = M/m, Vn = Tn\/m^. Eq. (IA.1|1 then 
becomes 

r1 _ _ _ _ 1 _ .X 

(A.36) 


Zn +-\<^+le 


Xn = 


The solutions thus only depend on c. The boundary conditions are as in Eqs. (IA.28f) . 
()A.29jl . (IA.30f) . for the symmetric and antisymmetric cases, respectively. 

This time we have only solved the eigenvalue problem numerically. The prohle 
A(z), the fundamental wave function i/olz), as well as the hrst excited symmetric wave 
functions, are shown in Fig. El for c = M/m = 1.4. Examples of eigenvalues are given 
in Table El 
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n 

smooth prohle, 

R = 2.5 
z/W a^/Jn{0) 

c = 1.4 

R = oo 

smooth prohle, 

.R = 3.15 
z/W a^/Jn{0) 

c = 2.0 

R = oo 

I'i"’ 

c 

0.000 

0.054 

— 

0.000 

0.034 

— 

0 

0.484 

0.134 

0.388 

0.033 

0.112 

0.030 

1 

3.020 

0.000 

3.016 

3.310 

0.000 

3.310 

2 

6.154 

0.112 

5.862 

5.517 

0.065 

5.513 

3 

9.430 

0.000 

9.373 

7.912 

0.000 

7.912 

4 

14.48 

0.112 

13.50 

11.16 

0.075 

11.15 

5 

18.51 

0.000 

18.17 

14.95 

0.000 

14.95 

6 

25.83 

0.112 

23.32 

19.27 

0.075 

19.21 

7 

30.17 

0.000 

28.92 

23.92 

0.000 

23.92 

8 

40.20 

0.112 

34.95 

29.20 

0.074 

29.03 

9 

44.63 

0.000 

41.37 

34.57 

0.000 

34.54 


Table 2: Eigenvalues for the smooth prohle. For the wave functions at origin one 
needs also the values of am, Ao/a, entering as a^jJn{0) oc (am)^/^(Ao/a)“^/^, if c, R are 
kept hxed; we have assumed am = 5.0/18.0, Aq/o = 8.75, for c = 1.4; and am = 
0.45, Ao/a = 25.0, for c = 2.0. 

Appendix B. Gauge invariant correlators in the Abelian case 

In the Abelian theory, Eq. (I2.1|l . the held strength tensor E/j, is gauge invariant. This 
allows one to measure directly various gauge invariant correlation functions displaying 
the essential features of the spectrum {m„}, as we will show with a specihc example. 
These correlators are however not available in the non-Abelian case. Other possibilities 
exist, but they contain either composite operators, making a qualitative distinction 
between genuinely conhning and Higgs-like phases difficult, or non-local operators, 
making the analysis of ultraviolet divergences as well as practical measurements hard. 
Let us dehne 

Op{x;z) = i (B.l) 

Ojr\x) = I^A{z)i,„{zWx;z) . (B.2) 

We might then consider, e.g., a weighted average over the z-direction, 

^ A( 2 :) JdxOp{x-,z)[Op{x+ r;z)]*, (B.3) 
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or, alternatively, a correlator only getting a contribution from some specified mode. 


Gf\r) = jdxO^^\x)[0^^\x + r)]\ 


(B.4) 


Employing the unitary gauge propagator (p, x, x' are 2d vectors) 



(B.6) 


one easily finds 




(B.6) 


n 



(B.7) 


Therefore, the spectrum {m„} again manifests itself in the form of the exponential 
decay. Note that in contrast to the force F{r;z), however, the wave functions 'ipni.z) 
do not appear in these predictions. 

Appendix C. The SU(Ac) string tension in 2d 

In this section we recall briefly the results for the static potential of pure SU(iVc) gauge 
theory in two dimensions. For a more detailed discussion see, e.g., ref. na. 

The result for the Abelian case is shown in Eq. (Q, with the 

non-Abelian case, that result is at leading order simply multiplied by an additional 
factor Ca = — l)/(2A^c), coming from the sum over the Hermitian generators of 

the fundamental representation, If we dehne Ca = 1 for U(l), the leading 

order result is then as written in Eqs. dSSI), (EED- 

It remains to show that there are no higher order corrections in the non-Abelian 
case. The naive argument goes as follows. As the potential is gauge hxing independent 
by construction, we can choose the gauge At = 0. But then the cubic and quartic 
interactions vanish, so that indeed no further corrections should arise. 
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